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1.0  Introduction 


The  overall  objective  of  the  program  is  the  photometric  and 
geometric  processing  of  UV  satellite  images.  The  data  source  is 
the  AIRS  scanning  UV  spectrometer  aboard  the  POLAR  BEAR  satel¬ 
lite.  The  processing  is  designed  to  provide  geographic-mapped 
images  which  describe  the  spatial  and  temporal  variability  of 
airglow  and  aurora  of  emitting  species. 

An  understanding  of  this  variability,  modeled  with  appro¬ 
priate  ground-based  measurements,  may  provide  maps  of  electron 
density  and  ionospheric  propagation  paths. 

The  AIRS  data  stream  provided  simultaneous  records  from  3 
detectors  received  at  several  ground  stations  from  several  passes 
per  day.  These  total  about  1500  useful  images  from  launch  to  Feb 
1987  when  the  satellite  lost  attitude  stability. 

The  objective  of  the  current  3-month  project  has  been  to: 

(  1  )  Evaluate  the  geometric  and  photometric  quality  of  the  AIRS 
data  collected  in  the  period  1  Jan  1987  through  28  Feb  1987. 

( 2  )  Provide  a  standard  processing  from  the  time-spaced  data 
stream  to  geographic  coordinates  of  a  photometrically-stan¬ 
dardized  image  which  will  clearly  show  the  relative  auroral 
activity  in  each  detector. 

13)  Provide  an  analysis  of  the  dependence  of  the  observed  count 
on  path-length  and  on  solar  flux  to  aid  in  modelling  the 
airglow  kinetics. 

The  major  progress  during  the  report  period  has  been 
development  of  programs  to  accomplish  the  following: 

(  1  )  Statistical  filtering  of  low-count  images  to  provide 
meaningful  brightness  estimates. 
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(2)  Use  of  observed  limb  locations  to  provide  attitude  correc¬ 
tion  in  daylit  images  and  attitude  estimation  in  night 
images  to  provide  accurate  ground  positions. 

(  3 )  Normalization  of  airglow  based  on  column  path  length  and 
solar  flux  which  allows  separation  of  aurora  from  airglow 
portions  of  the  images. 

(4)  Multi-band  presentation  of  auroral  features  which  clearly 
define  differences  between  the  location  of  emitting  species. 

{ 5 )  Batch  processes  which  will  take  a  given  record  from  raw 
input  to  final  mapped  aurora  and  airglow,  and  analysis  of 
spatial  dependence. 

This  report  presents  a  brief  summary  of  each  topic  followed 
by  a  series  of  Appendices  which  provide  more  detailed  discussion 
of  each  topic. 

Appendix  1  presents  a  summary  of  programs  developed  during 
the  report  period. 
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2.0  Data  Sources 

The  image  data  is  merged  with  the  orbital  position  and  atti¬ 
tude  data  into  a  specified  format  on  a  PC  diskette,  which  is  the 
raw  input  for  the  processing  system  described  below. 

The  raw  image  data  in  each  detector  is  a  count  of  sub-satel¬ 
lite  column  brightness  sampled  at  7  millisec  intervals  during  a 
135  degree  mirror  scan  around  nadir  over  a  3  second  period.  Each 
detector  isolates  emission  in  a  narrow  wavelength  band  character¬ 
istic  of  a  distinct  emitting  species  and  process. 

This  scan  presents  a  horizon-to-hor izon  record  of  4000  km 
across  the  orbital  path.  A  typical  record  from  a  station  during 
a  satellite  pass  includes  about  200  scan  lines  taken  over  600 
seconds  representing  a  4000  km  along-orbit  track. 
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The  program  RAWCVT  has  been  rewritten  to  convert  this  raw 
diskette  into  a  240  line  x  240  column  image,  an  orbital  data  file 
and  a  wavelength  file.  The  program  copes  with  data  gaps  and 
tests  whether  the  raw  data  is  packed.  If  packed,  program  UNPAK 
is  invoked  to  generate  a  non-saturated  image. 
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Appendix  2  lists  details  of  this  program. 

3.0  Statistical  filtering  of  low-count  images 

Count  rates  in  all  detectors  were  usually  in  the  0  to  10 
range  and  rarely  exceeded  20.  This  data  is  best  treated  as  a 
Poisson  process.  A  count  of  1  has  a  100  pc  error  and  a  count  of 
10  a  30  pc  error.  The  error  can  be  reduced  to  10  pc  by  'lumping' 
samples  until  100  counts  are  accumulated,  effectively  trading 
spatial  resolution  for  a  statistically-acceptable  error.  A  pro¬ 
gram  FILTRAW  was  developed  which  provides  the  necessary  tradeoff 
at  each  point  of  the  image,  while  preserving  the  original  per¬ 
centile  distribution. 

Appendix  3  lists  details  of  this  program. 

4.0  Attitude  data 

During  observation  when  the  satellite  is  in  sunlight,  atti¬ 
tude  is  derived  from  a  magnetometer  and  solar-aspect  sensor  with 
about  +  2  degree  accuracy  and  reported  for  each  3-second  scan. 
Pitch,  yaw  and  roll  showed  5-15  minute  cycle  amplitudes  ranging 
from  0  to  several  degrees. 

At  other  times,  no  attitude  data  is  available  from  the 
satel 1 ite . 


However,  a  bright  limb  is  observed  from  the  emitting  layer 
and  a  program  ROLADJ  has  been  developed  to  estimate  the  roll 
required  to  place  the  limb  at  a  defined  altitude.  In  several 
images,  both  limbs  are  visible  and  defined  the  emitting  layer 
height  as  1 50  +  20  km. 


The  roll  defined  by  the  limb  showed  about  a  1.5  degree  off¬ 
set  from  the  available  magnetometer-reported  rolls. 

The  average  roll  adjustment  for  each  of  the  cases  tested 


Code 


10  11  12  13  14 


Rol-adj  1.7  1.0  1.8  2.0  2.2  1.5  1.6  2.8  3.6  1.0  2.3  5.0 

Appendix  4  lists  details  of  this  program. 

5.0  Geometric  mapping  of  the  raw  image 

The  observed  count  came  from  the  line-of-sight  determined  by 
mirror  scan  angle  and  satellite  attitude.  At  angles  near  nadir, 
this  is  a  near-vertical  column  from  satellite  to  ground.  As  the 
scan  angle  increases,  the  column  is  more  oblique  and  above  60 
degrees  from  nadir  does  not  reach  the  ground  but  continues  into 
space.  Near  the  left  and  right  edges  of  the  raw  image  where  this 
condition  applies,  some  samples  are  not  intercepting  the  ground 
and  cannot  be  mapped  onto  earth  coordinates. 

At  angles  where  the  line-of-sight  scan  intercepts  the 
ground,  there  is  a  single  oblique  intercept  with  the  emitting 
layer  at  about  150  km,  which  will  be  closer  to  satellite  nadir 
than  the  ground  intercept  distance. 

Program  GEOMIM  and  SOLEL  have  been  written  to  compute  for 
each  image  point: 

( 1  )  The  geographic  coordinates 

(2)  The  local  solar  flux  =  sin  (solar  elevation) 

( 3 )  The  oblique  path  length  relative  to  vertical  which  reaches 
about  5  X  the  vertical  path  at  the  ground  tangent.  Raw 
samples  at  larger  scan  angles  cannot  be  assigned  ground 
coordinates . 

A  program  WARPIM  has  been  written  to  'warp'  the  image  from 
the  input  sample  coordinates  to  ground  coordinates. 
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Appendix  5  lists  details  of  this  program. 

6.0  Normalization  of  Airqlow 

A  program  GEOMLSQ  has  been  written  to  test  the  dependence  of 
airglow  intensity  on  the  path  length  and  solar  flux. 

The  program  found  the  dependence  of  count  on  some  power  of 
path  and  some  power  of  solar  flux. 


cnt  =  CO  *  (path  **  -Cp)  *  (sin(solel)  **  Cs ) 
in  sunlight  by  a  least  square  fit  of 

loglO(cnt)  =  CO  +  Cp  *  log  1 0  (  1 /path )  +  Cs  *  log10(sin  sol 
elev ) 

and  in  dark  by 

loglO(cnt)  =  CO  +  Cp  *  log  1 0 ( 1 /path ) 

Detector  2  (1356  A)  was  the  only  group  with  a  count  which 

did  not  require  very  heavy  smoothing.  In  this  group,  the 
observed  count  of  non-auroral  sunlit  points  was  well-fitted  to  an 
expression  proportional  to  path  length  ( opticAlly-thin )  and  to 
sqrt  (solar  flux).  The  non-sunlit  airglow  did  not  show  this  path 
dependence  ( optically-thick? ) .  The  other  detectors  were  so  weak 
that  no  r<  liable  fit  was  possible.  These  facts  should  be  useful 
in  defining  the  kinetics  of  the  emission  process. 


Table  1  shows  the  outputs  of  GEOMLSQ  on  the  cases  listed  for 
each  detector  for  sunlit  and  non-sunlit  portions  of  the  image. 

Table  column  headings  are: 
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Table  1  Summary  o-f  geomlsq  results  on  13  scenes 
prepared  by  REPLSQ 
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0 

28 

22.3 

-0.  17 

0.00 

0 .  1 30 

94 

C' 

539 

12 

024 

23.62 

0 

31 

21.6 

0.35 

0.00 

0.  220 

91 

C> 

977 

13 

029 

23.  27 

0 

39 

10.  9 

-0.06 

0.00 

0.  260 

lOO 

1026 

14 

032 

23.05 

0 

34 

12.5 

-0.51 

0 . 00 

0,  320 

91 

C> 

□73 

15 

04  3 

22.  30 

0 

33 

26.5 

0.82 

0.00 

0.  270 

69 

0 

3023 

kv 


COi  count  in  amplified  filtered  image  normalized  to 

path  length  1  and  in  sunlit  cases  to  a  sin(solar  elev)  of 

0.01 

Cp  Isq  fit  i/path  dependence 

Cs  Isq  sin( sol_elev  )  dependence 

sgO  stdvn  of  input  points  in  loglO  units  (0.13  equiv  to 

+  30%  ) 

pcvr  percent  of  variance  remaining  after  fit  (i.e.  10  pc 

equivalent  to  reduction  of  stdvn  to  30  pc  of  input) 

Note  that  if  the  pcvr  is  above  50  pc,  the  data  scatter  is  so  high 
that  there  is  no  significant  dependence  on  the  parameters 
pcvx  percent  of  variance  remaining  after  fit  forcing  Cp  to  -1 
and  Cs  to  0.5 

Note  that  if  pcvx  is  close  to  pcvr,  the  force  fit  is  as  valid  as 
the  ' f ree '  Isq  fit 

Ray  the  Rayleigh  count  normalized  from  COl  correcting  for 

amplification  from  raw  to  filtered  count  and  the  detector 
sensitivity  ( Ray leighs/orig  count) 

Appendix  6  lists  details  of  this  program. 

7.0  Multi-band  presentation  of  auroral  features 

As  noted  above,  it  is  possible  to  define  the  airglow 
relationship  to  solar  flux  and  path  length.  However  auroral 
features  do  not  match  the  expressions  and  may  be  'feature- 
extracted'  from  the  airglow  and  separately  presented. 

The  relative  intensity  of  local  aurora  in  two  bands  may  most 
clearly  be  shown  after  a  histogram-equalization  program  HISEQ  is 
performed.  By  this  process  the  lowest  1  percentile  is  converted 
to  relative  intensity  1,  the  50  percentile  to  50,  and  the  99 
percentile  to  99.  In  practice  the  relative  intensities  are  scaled 
0-255  . 
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The  two  bands  can  be  co-displayed  by  a  program  MULBND  as  a 
2-color  image  with  one  band  in  red  and  the  other  in  green.  If  the 
two  bands  match,  a  yellow  image  of  a  given  intensity  will  be  dis¬ 
played,  else  a  tendency  to  red  or  green  will  be  displayed.  In 
practice,  the  program  adds  a  'blue'  signal  equal  to  the  weaker  of 
the  two  signals  so  the  matched  color  'yellow'  is  seen  as  shades 
of  grey  from  black  to  white.  (The  problem  of  subtracting  one 
detector  signal  from  another  as  a  'contaminant'  has  not  yet  been 
addressed . ) 

Appendix  7  lists  details  of  this  program. 

8 . 0  Batch  processes 

Mulp.bat  carries  a  given  detector  through  all  steps  required 
by  the  flowsheet  after  roladj  to  produce  labelled  images  as  disk 
files.  Preparation  of  a  full  group  of  32  display  images  from  a 
pair  of  input  detectors  takes  about  15  minutes.  Required  inputs 
for  Mulp.bat  are  a  pair  of  raw  image  files  and  data  files 
rawNN.d,  wvlNND.d,  and  rajNN.d  roll  adjustment. 


Muld.bat  displays  the  labelled  images  for  demonstration  pur¬ 
poses.  After  photography  of  appropriate  images,  the  image  files 
should  be  deleted  because  of  heavy  disk  usage. 

Appendix  8  lists  details  of  these  batch  files. 

9.0  Photographs 


Photographs  from  the  monitor  screen  demonstrating  various 
stages  in  the  processing  are  reported  in  Appendix  9. 
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APPENDIX  1.  Programs,  Data  base  structure,  and  Flowsheets 
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Appendix 

1  -a . 

Program 

Type 

BEEP 

Util 

CB 

Util 

CBDIAG 

Util 

CLEANIM 

AIR 

DIFFIM 

Img 

DIRFMT 

Util 

DIROS 

Util 

DIRDUP 

Util 

FFTAUR 

AIR 

FILTIM 

Img 

FILTRAW 

AIR 

FINDX 

Util 

FULLDIR 

Util 

GEOMIM 

AIR 

GEOMLSQ 

AIR 

GRAF 

Util 

GRAFIM 

Img 

GRIDIM 

Img 

HISEQ 

Img 

INITIM 

Img 

LANDMASS 

Img 

MEDIAN 

Img 

MOUSE 

Img 

MODIM 

Img 

MOVEIM 

Img 

MULBND 

AIR 

NEWSCL 

Img 

OVERLIM 

Img 

PATCH 

Img 

PAUSEX 

Util 

POISCMP 

AIR 

POISEST 

AIR 

POISSON 

Img 

PRINTIM 

Img 

RAWCVT 

AIR 

RAWFIT 

AIR 

REFMT 

Util 

REPAK 

Img 

REPLSQ 

AIR 

ROLADJ 

AIR 

ROTIM 

Img 

SNAPIM 

Img 

SNAPINT 

Img 

SOLEL 

AIR 

TEXTHDR 

AIR 

TEXTIM 

Img 

UNPAK 

AIR 

WARPIM 

AIR 

ZOOMIM 

Img 

New  or  modified  source  programs 


Beeps  to  console 
Beautify  C  program 
Diagnose  C  program 

Remove  isolated  high  noise  points  from  image 
Difference  between  2  images 
Modify  dir  listing 

Modify  dir  listing  to  provide  bat  file  for  del  or 
copy 

Find  dup  files  in  full  dir  listing 
Test  for  amplitude  in  aurora  after  filter 
Filter  smoothing  with  possible  threshold 
Narrow  and  wide  smooth  combine  to  filtered  image 
Find  and  print  lines  not  containing  phrase 
Generate  full  listing  of  disk  from  'TREE'  listing 
Generate  matrix  of  ground  coord  for  raw  image  points 
Least  sq  fit  of  brightness  to  path  len  and  sol  flux 
Graph  input  data  set 

Graph  to  image  selected  image  line  brightness 
Overlay  grid  on  image 
Histogram  equalize  image 
Initiate  image  board 

Overlay  landmasses  and  lat/lon  cord  on  image 
Median  image 

Provide  print  around  moveable  cursor  mouse  on  image 
Stretch  rotate  1 -d  smooth  image 

Move  image  from  qd  to  qd  or  disk  to  qd  or  qd  to  disk 

Two-color  output  from  2  input  bands 

Rescale  image  brightness 

Overlay  image  ( landmass )  onto  image 

Patch  noisy  segments  from  image 

Pause  in  program 

Compute  Poisson  statistics 

Estimate  Poisson  statistics 

Generate  Poisson  statistics 

Print  vals  from  image 

Convert  Raw  diskette  to  image  file 

Reprocess  raw  orbital  parameters 

Reformat  file 

Pack/unpack  log  expansion 

Summarize  all  Isq  files 

Compute  correct  roll  from  limb  location 

Rotate  or  squish  image 

Snaps  to  quadrant 

Snaps  and  integrates  to  quadrant 

Compute  solar  elev  for  each  point  in  image 

Extract  text  from  hdr  file  and  overlay  for  demo 

Label  Image 

Unpack  packed  satellite  data 

Warp  raw  image  to  ground  coordinates 

Zoom  image  to  quad  or  1/4  quad 


Appendix  1-b.  Scene  Processing  Flowsheet  for  AIRS 


The  flowsheet  below  shows  processes  and  input  and  output  files 
generated  for  a  given  scene.  A  scene  NN  (coded  by  a  2-character 
symbol,  01 -ZZ  possible)  is  defined  as  a  multi-detector  record  for  a 
given  site,  day  and  pass. 

The  geographic  data  files  are  required  for  processing  each 
individual  detector  image.  These  data  are  common  to  all  detectors  and 
are  generated  by  this  flowsheet. 


INPUTS 

PROCESS 

OUTPUTS 

Data  Images 

Data 

Images 

RAW 

1  rawcvt 

1  RAWNN.d 

RAWNND. i 

DISKETTE 

1 

1 

!  WVLNND.d 

RAWNN.d 


FILNN.d 

HDRNN.d 


RAWNN.d 

RAJNN.d 


FILNN.d 


GEONN.d 


GEONN.d 

SELNN.d 


HDRNN.d 


!  rawfit 

!  FILNN.d 

1 

i 

1 

!  HDRNN.d 

1  roladj 

1 

1 

I  RAJNN.d 

1 

1 

1 

1  rawfit 

!  FILNN.d 

1 

1 

1 

!  HDRNN.d 

1  geomim 

1 

1 

!  GEONN.d 

1 

t 

1 

1 

1  solel 

I 

!  SELNN.d 

1 

1 

t 

1 

I  warpim 

1 

1 

1 

1 

I 

SELNN. i 

1 

f 

1  landmass 

1 

1 

LHNN . i 

i  10 

✓ 
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Appendix  1 -c .  Detector  Processing  Flowsheet  for  AIRS 

The  flowsheet  below  shows  processes  and  input  and  output  files 
generated  for  photometric  processing  for  each  useable  detector  in  a 
given  scene.  The  individual  detector  files  are  coded  NND  where  D  is 
expected  to  have  a  value  1  to  4 . 

The  geographic  data  files  are  combined  with  the  detector  files  to 
provide  geographically-oriented  5000  x  5000  km  scenes. 

INPUTS  PROCESS  OUTPUTS 

Data  Images  Data  Images 

RAWNND.i  I  unpak  !  UNPNND.d  UNPNND.i 


UNPNND.i 

1  filtraw 

1 

I 

1 

1 

1 

1 

FILNND. i 

GEO  NN.d 
GEO  NN.d 

FILNND. i 

1 

1 

1  warpim  opi 

1 

1 

!  AMPNND.d 

1 

1 

HDRNND.d 

SELNN.  i 

1 

!  geomlsq 

!  LSQNND.d 

AGLNND. i 

AMPNNd.d 

UNPNND.d 

WARNND.  i 

1 

1 

( 

1 

1 

1 

AURNND. i 

AURNND. i 

!  hiseq 

1 

1 

1 

t 

1 

1 

AUHNN.i 

AUHNN1 . i 
AUHNN2.i 

LHNN.  i 

1 

1 

1  mulbnd 

1 

1 

i 

f 

1 

1 

1 

1 

AURNN 1 2 . i 

AUHNN1 . i 

!  overlim 

1 

1 

DISNNl . i 

AUHNN2 . i 

1  texthdr 

1 

1 

DISNN2 . i 

AURNN1 2 . i 

— 

DISNNl 2 . i 

photographs 
for  publication 


any  image 

1  grafim 

1  graphed  count  along 

1 

I  sampled  lines 

photographs  for  publication 
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APPENDIX  2. 


Conversion  of  raw  diskette  to  image  and  data  files 


RAWX . BAT  file  to  process  several  groups  of  input  diskettes 
goto  aa 

del  \air\raw1 3 .d 

pausex  insert  diskette  sn87029.d3  for  raw133 
rawcvt  s87029  133  es4420  <  a:sn87029.d3 
moveim  15  \air \raw1 33 . i 

pausex  insert  diskette  sn87029.d2  for  raw132 
rawcvt  S87029  132  es4420  <  a:sn87029.d2 
moveim  15  \air \raw1 32 . i 

pausex  insert  diskette  sn87029.d1  for  raw131 
rawcvt  S87029  131  es4420  <  a:sn87029.d1 
moveim  15  \air \raw 1 3 1 . i 
:  aa 

del  \air\raw14.d 

pausex  insert  diskette  sn87032.d3  for  raw143 
rawcvt  S87032  143  es16050  <  a:sn87032.d3 

moveim  15  \air \raw1 4 3 . i 

pausex  insert  diskette  sn87032.d2  for  raw142 
rawcvt  S87032  142  es16050  <  a:sn87032.d2 

moveim  15  \air \raw 1 42 . i 

pausex  insert  diskette  sn87032.d1  for  raw141 
rawcvt  S87032  141  es16050  <  a:sn87032.d1 
moveim  15  \air\raw1 4 1 . i 

RAWCVT  help 

USAGE  rawcvt  SYYDDD  NND  ( raW  bsB  esF  ssS  fsS  IpL  hpH )  <  ASCII 
EXAM  rawcvt  S86338  Oil  <  A : S87008h2 . dat 
INPUTS  ASCII  file  with  original  BC  file  format 
OUTPUTS  in  \AIR  directory 

rawNN.d  if  not  exist  else  check  vs  rawNN.d 
wvlNND.d  with  wvln  of  detector  D  image 
qd=2  xfered  to  qd=l  (rot  if  necessary) 

SYYDDD  converted  in  rawfit  to  site  yr  and  day 
NND  AIR  code  and  detector 

stretches  input  to  image  (240  line  x  240  col) 
estimates  lines  from  file  header  start  and  final  sec 
if  conflict  report  start  and  final  sec  in  line  hdr 
can  rerun  with  options  ssS  and  fsS  described  below 
save  image  by  'moveim  15  rawNND.i' 
options  if  present 

raW  W= 1  output  q1  raw  data  326  cols  x  input  lines 
bsB  flags  begin  sec  for  image  (if  bad  section) 
esE  flags  end  sec  for  image  (if  bad  section) 
ssS  S=  start  sec  to  override  header  start  sec  error 
fsS  S=  final  sec  to  override  header  final  sec  error 
IpL  hpH  low  high  input  kol  with  non-zero  (def=5  330) 
bgB  B=2  quits  after  hdr  read 


Terminal  print  from  RAWX.BAT 

C : \B0B1 >pausex  insert  diskette  sn87043.d3  for  raw153 
insert  diskette  sn87043.d3  for  raw153  continue.. 

C : \BOB1 >rawcvt  s87043  153  es4l90  <  a:sn87043.d3 
C:\B0B1\BIN\RAWCVT.EXE  s87043  161  es4190 
opening  junkraw  for  diagnostics 
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opening  new  scene  file  \air\rawl6.d 
opening  detector  file  \air\wvl1 61 -d 
writing  wvl  1 596 

raw  SS=3805.1  FS=4501.1  sec=696.0  lines=233 
req  BS=3805.1  ES=4190.0  sec=384.9  lines=129 
adjusted  ESEC  to=  4189.10 

final  write  at  CSEC=4189.1  lin=128  lwp=238  lwr=240 

ilin  129  Iwr  240  ESEC=4189.1  CSEC-ESEC=  0.0 

report  min(9,cnt)  for  val=  0-60  if  unpacked  33  35  37  zero 

val  0  5  10  15  20  25  30  35  40  45  50  55  60 

cnt  9999999999993322043000101200000000000000000000000000000000000 

No  errors  in  rawcvt 

Rotate  image  for  N  up 


File  junkraw  generated  by  rawcvt 


raw  SS=3805.1  FS=4501.1  sec=696.0  lines=233 
req  BS=3805.1  ES=4190.0  sec=384.9  lines=129 
adjusted  ESEC  to=  4189.10 


raw 

out 

rawNN 

curr 

imag 

read 

true 

calc 

lin 

Iwr 

prnt 

CSEC 

RSEC 

dsec 

1  rd 

1  wr 

0 

0 

1 

3805  .  1 

0.0 

0.0 

0.0 

0.0 

1 

2 

1 

3808. 1 

3.0 

3.0 

1  .0 

1  .9 

2 

4 

1 

3811.1 

6.0 

3.0 

2 . 0 

3.7 

5 

9 

2 

3820 . 1 

15.0 

3.0 

5 . 0 

9.3 

125 

233 

38 

4180.1 

375.0 

3.0 

125.0 

232.6 

1  26 

234 

38 

4183.1 

378.0 

3.0 

126.0 

234.4 

1  27 

236 

38 

4186.1 

381.0 

3.0 

127.0 

236.3 

1  28 

238 

39 

4189.1 

384.0 

3.0 

128.0 

238.2 

final 

write  at  CSEC=4189.1  lin 

=  128 

lwp=238  lwr=240 

ilin 

1  29 

Iwr  240 

ESEC=4189.1  CSEC-ESEC=  0.0 

report  min(9,cnt) 

for  val=  0-60 

if  unpacked  33 

1  35  37 

val  0 

5  1 0 

1  5 

20  25 

30 

35  40 

45 

cnt  9999999999993322043000101200000000000000000000000000000000000 
No  errors  in  rawcvt 
Rotate  image  for  N  up 


APPENDIX  3.  Amplif ication  and  filtering  of  raw  image 

USAGE  filtraw  ( nW  wW  sS  pP  f F ) 
inpuv  ql  amplifies  and  smooths 
outputs  q2=aurora  q3=airglo  q4=final 
reports  final  amplification  in  file  junkamp 

pP  P=pcaur  pet  image  with  narrow  smth  wid  (def=25) 
nW  W=narrow  smoothing  width  (for  aurora) 
wW  W=wide  smoothing  width  (for  airglow) 
sS  number  smths 

fF  F^final  value  wanted  at  pcaur 
defaults  n3  wl9  si  fSO 
if  set  nO  wO  sO  fO  selects  n  w  s  f 
EXAM  filtraw 

EXAM  filtraw  n3  w9  si  p)0 

Comparison  of  image  count  with  Poisson  distribution 

A  Poisson  distribution  is  expected  for  a  random  process 
where  the  sample  rate  is  high  compared  to  the  count  rate. 

The  expression  for  probability  of  count  N  in  a  Poisson 
process  is 

P(N)=  SAV  **  N  /  N!  *  exp(-SAV) 

with  both  standard  deviation  and  signal/noise  ratio  (SNR)  = 
sqrt(SAV  ) 

As  an  example,  if  a  count  rate  of  100/sec  is  sampled  at  a 
rate  of  100/sec  the  average  count  per  sample  (SAV  =  tot 
counts/tot  samples)  will  be  1  with  a  distribution  of  about  37% 
with  count  0,  37%  with  count  1  and  18%  with  count  2.  (The 
Poisson  expression  approaches  a  Gaussian  distribution  at  high  SAV 
values.  )  Table  1  presents  selected  values  of  Poisson  statistics. 

Table  i.  Poisson  probabilities  for  selected  sample  average  count. 

SAV  sample  average  count 

P(N)  1000  *  probability  of  count  N 


1  1  5 


;o  2  5 


0.01  P(N)=  990  9 

0.03  P(N)=  968  30 

0.1  P( N  )  =  904  90  4 

0.3  P(N  )=  728  230  36  3 

1  P(N)=  367  367  183  61  15 
3  P(N)=  42  133  211  223  176  1 


10  P(N  )=  0 

20  P(N)=  0 


58  26  10 


18  37  63  90112125125113  34 


10  51  88  34 


f.'. 


To  achieve  reasonable  SNR  =  sqrt(SAV)  it  is  possible  to  tradeoff 
time  (or  spatial)  resolution  by  increasing  the  sample  period,  or 
equivalently  to  sum  enough  samples  to  increase  the  sample  average 
SAV.  Specifically,  the  SNR  =  3  if  enough  samples  are  summed  to 
bring  SAV  to  10  and  SNR  =  10  if  SAV  is  brought  to  100. 

AIRS  images 

Are  the  low-count  AIRS  images  typical  Poisson  processes? 
Table  2  presents  the  statistics  for  a  night  AIRS  image  in  each  of 
3  bands.  The  data  confirm  that  at  the  low-count  region  of  the 
images  are  closely  approximated  by  Poisson  distributions. 

This  implies  that  variable-area  summing  must  be  performed 
using  different  summing  areas  for  different  SAV  regions.  The  area 
will  be  100/SAV  pixels  for  SNR  =  10.  Such  a  summing  is  in  fact  a 
smoothing  (or  averaging)  process  and  each  summing  area  leads  to  a 
different  image. 

In  practice,  two  smoothing  areas  are  sufficient  for  most 
purposes.  A  large-smoothing  area  image  based  on  the  low-count 
region  and  a  smaller-smoothing  area  image  based  on  the  high-count 
area.  The  final  output  uses  the  value  from  the  narrow-smoothed 
image  if  it  is  above  the  narrow-smooth  threshold,  else  from  the 
wide-smoothed  image. 

Notes  on  the  smoothing  process 

The  original  image  is  first  'amplified'  by  a  constant 
determined  from  the  input  histogram  to  bring  the  99  percentile 
value  to  a  value  of  200.  This  is  necessary  because  averaging  a  50 
pixel  area  containing  90  percent  O's  and  10  percent  1's  will 
yield  zero  but  if  the  1's  have  been  amplified  to  50's,  a 
meaningful  value  of  5  will  be  obtained. 

It  is  more  efficient  to  carry  out  two  1 -dimensiona 1 
smoothings  than  an  area  smoothing.  A  horizontal  (  line-by-11 ne > 
smoothing  of  input  is  followed  by  a  vertical  smoothing  of  the 
first  output.  The  smoothing  width  in  each  direction  will  be  the 
sqrt  of  the  desired  smoothing  area. 
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Table  2 


Comparison  of  Poisson  and  AIRS  Scene  Count  distribution  in  non- 
auroral  regions  of  a  nightime  scene. 
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SAV 

= 

2  . 

63  N 

50 

=  2. 

90 

p 

I  0* 

n 

0 

23 

29 

33 

4  1 

47 

52 

57 

62 

66 

71 

75 

79 

s 

1  0* 

n 

0 

22 

27 

32 

40 

48 

56 

65 

76 

91 

1  1  2 

1  38 

171 

'T'  b 

/  Ex 

- 

95 

93 

96 

97 

1  02 

1  07 

1  1  4 

1  22 

1  37 

1  57 

184 

216 

S xplanation  of  symbols: 


is  the  cumulative  percent  with  count  >=  N. 
to  low  values  at  high  N 


Ranges  from  1000  at  N=0 


Iqf  r 

w  1  n 
?  ■' 

N  5  0 


log  fraction  remaining  reported  as  100  *  loglO  (100/cmpc).  N  is  a 
linear  function  of  Igfr  for  small  SAV 

Wavelength  of  image 

Sample  average  =  total  counts/total  samples 
Interpolated  N  for  50  cumulative  percent  of  samples 


P_10*n  The  expected  Poisson  N  interpolated  to  selected  Igfr  and  reported 
as  10  *  N . 

S  i0*n  The  scene  distribution  of  N  interpolated  to  selected  Igfr  and 
reported  as  10  *  N. 

ijb, Ex  The  ratio  of  Observed  Scene  to  Expected  Poisson  N  reported  as 

percent  Scene  N/Poisson  N 


! j 'it pur  from  filtraw 


':pt  ions  - 

nsm  number  of  smooths  in  each  direction 

li^wid  narrow  smoothing  width 

w  wid  wide  smoothing  width 

pcau  percent  aurora 

(pc  image  with  narrow  smoothing  width) 
lafr  au  log  cum  fraction  at  pcau  as  -  1 00 . *log1 0(pcau/ 1 00 .  ) 

50pc  equiv  30  25pc  equiv  60  1 Opc  equiv  100  Ipc  equiv  200  fin_aur 
final  value  in  output  image  for  pcaur 

16 


r 


options  nsm=l  n_wid=3  w_wid=19  pcau=25  lgfr_au=  60  fin_aur=50 


histog  report- 

input  avg=  6.892  avg/av0=  1.00  highest  val=56  iaur=9 
input  avg=6.892  aur  value  =9.00  mulau=5.6  to  bring  fin_aur  to  50 

lookup  table  out  for  inp 


inp 

0 

1 

2 

3 

4 

9 

40 

45 

80 

1  20 

1  60 

200 

240 

lut 

0 

6 

1  1 

1  7 

22 

50 

222 

250 

255 

255 

255 

255 

255 

for  1 

given 

cumulative 

per  1 

Tii  1  . 

above 

_ 

cmpm 

1  000 

630 

501 

398 

316 

251 

1  58 

1  00 

63 

39 

25 

1  5 

1  0 

6 

3 

Igf  r 

0 

20 

30 

40 

50 

60 

80 

1  00 

1  20 

140 

1  60 

180 

200 

220 

240 

inp 

0 

2.7 

4 . 7 

6  .  1 

7.6 

9 

1  2 

1  7 

24 

32 

36 

42 

46 

46 

46 

v*rnu 

0 

1  5 

26 

34 

42 

50 

69 

95 

1  34 

1  76 

200 

232 

255 

255 

255 

average  of  values  above  pcau  and  estimated  smoothing  width  for 
stdvn  =10  pc 

avau=  20.067  est  n_wid=2  using  n_wid=3 

average  of  values  below  pcau  and  estimated  smoothing  width  for 
stdvn= 1 0  pc 

avag=  3.284  est  w_wid=7  using  w_wid=19 


suggested  nsm=l  using  nsm= 1 
result  of  smooth  and  read justments- 

first  smth  199=237  avg=36.5  iaur=46  anor=1.09  mulau=6.04  iaur=49 
output  avg=  42.105  avg/av0=  6.11  highest  val=255  iaur=51 

ratio  of  out/ (  inp*a:np  ) 

val  1  24  30  36  42  52  69  99  143  182  209  239  253  255  255 

v/vO  -  160  115  105  100  104  100  104  106  103  104  103  99  100  100 

filtraw  amplification  is: 
amp_f ac=6 . 1 09  av_inp=6.892  av_out=42 . 1 05 

Note  the  aiuplif  ication  finally  used  is  based  on  unpacked  vs 
filtered  values  in  the  area  warped  onto  the  ground. 


Comparison  of  input  and  filtered  image 


input 


col  = 

89 

1  00 

1  1  0 

1  2  1 

131 

142 

1=101 

1 

0 

2 

0 

4 

9 

1  =  109 

3 

3 

1  3 

8 

6 

9 

*  117 

8 

1 

5 

3 

3 

5 

1=12  5 

4 

6 

4 

7 

7 

2 

1  --  1  3  3 

5 

1 

2 

3 

1 

3 

52 

1  filtered 

1  89100110 

1  2  1 

131 

1  42 

1  52 

1  6 

1  0 

1  3 

1  9 

26 

32 

68 

74 

4 

20 

24 

7  1 

38 

39 

57 

34 

2 

30 

32 

35 

35 

32 

28 

24 

2 

34 

31 

28 

26 

25 

24 

23 

3 

32 

29 

25 

24 

24 

24 

24 
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APPENDIX  4. 


Roll  Ad]ustment 


SAGE:  roladj  NN  (AIM  IIL  IhL  kwK  mxM  muM  ipR  nsM  prP  rrR  yrY  buG 
htH  nsM  dgD  +  possible  overrides  psP  ysY  rsR  ksK  tsT) 

EXAM  roladj  01 

EXAM  roladj  01  lllOO  lh200  mx 1 

EXAM  roladj  01  ?  (interactive  query) 

INPUTS  hdrNN.d  filNN.d  and  image  rawNND.i  in  ql 


\AIR\ra jNN.d 


Eyeball  raw  images  to  find  detector  with  best  limb  and  assign  11  Ih 
jjrora  area  to  exclude  a  helpl  for  full  help 

Extraction  of  Adjusted  Roll  from  Image 
The  roladj  program  (Limb  Horizon) 

The  program  has  the  objective  of  determining  an  adjusted 
roll  angle  which  will  place  an  observed  bright  limb  at  a  defined 
a  Ititude . 

D 3 0 a  Processing 

For  AIRS  satellite  total  scan  angle=130.4  deg  The 
deg/ image_col=  0.54  from  130.4  scan_deg/240  sample  columns.  From 
the  satellite_earth  geometry  for  satellite  alt=1000  km.  The  scan 
angle  from  nadir  to  the  satellite  tangent  to  ground=  60.0  deg. 
The  scan  angle  from  nadir  to  the  satellite  tangent  to  limb  as  a 
function  of  limb  height  and  the  distance  in  columns  from  limb  to 


ground 

tangent 

is  as  follows: 

Relation 

of  limb 

height  to  limb  limb 

w  :  d 

^  h  in 

degrees 

and  image  columns 

t  o  t 

deg 

col 

gnd 

119.6 

220.0 

Imc 

each_side 

h  U 

col  Imb  gnd 

'22.8 

226.0 

3.0  ” 

1  2  0 

12  3.4 

2  2  7.2 

3 . 6 

;  50 

124.4 

229.0 

4  .  5 

’  7  0 

12  5.1 

2  30.2 

5  .  1 

’  9  0 

125.8 

231.5 

5.8 

2  50 

127.9 

2  3  5.3 

7  .  7 

In  all  cases  where  2  or  more  wavelengths  showed  limbs,  the 
limbs  were  located  within  ^  1  column. 

In  all  cases  where  both  limbs  were  present  the  limb-limb 
distance  was  230  +_  1  column  so  a  'best  height'  estimate  is  170  + 
2  0  km . 


Terminal  Output  from  rolladj 

C:\BOB\BIN\ROLADJ.EXE  06  mxl 
hdr  0  AIR  06 

hdr  1  Site  SON 

hdr  4  LMT  9.86 

hdr  9  Radj  0.0 

opt ions- 

right  limb  height=150  kol_limb=239  11=240  lh=0  mu=2  ipr=1  nsm=3  dg=0 

smoothing  image  from  q2  to  q2  vertically 

found  strongest  limb  line=239  kol=234  britness=161 


scan  limb  limb 

ang  ht  col 

130.4  150  235.1  0.00 


lm_gd  deg/  sat_ion  sat_gnd  pit  yaw  roll 

cols  col  angle  angle 

4.38  0.54  61.93  59.55  0.9  0.6  -3.9 


Note  roll  is  positive  for  right  wing  down. 

calculation  of  roll  from  right  limb  location 
lin  p  y  roll  limb  ang  limb_col 

inp  out  adj  1ft  rgt  1ft  rgt  wid 


0 . 1 *km  gnd  I 
grnd_col  I  nadirl 
1ft  rgt  wid I  1ft  rgt! 


-2.2  1.6  57.2  61.7 


calc  found 
10.0  238.0  228.0 


233.6 


0  200 


1  8 

1 

3 

-3.8 

CM 

CM 

1 

1  .6 

57 . 2 

61.7 

10.0 

238.0 

228 . 0 

14.7 

233.6 

218.9 

200 

200 

36 

1 

8 

-3.4 

-1.3 

2  .  1 

58.  1 

60.8 

8.3 

236.3 

228 . 0 

13.0 

231.9 

218.9 

1  99 

1  99 

54 

1 

8 

-2 . 6 

-0.3 

2 . 3 

59.1 

59.8 

6.4 

234.4 

228 . 0 

11.2 

230.0 

218.8 

1  99 

1  99 

72 

1 

3 

-1.9 

1  .  0 

2.8 

60.4 

58.5 

4.0 

232.0 

228 . 0 

8.9 

227.6 

218.7 

200 

201 

90 

1 

-3 

-1.8 

0.4 

2 . 2 

59.9 

59.0 

4.9 

232.9 

228 . 0 

9.8 

228.5 

218.7 

201 

200 

1  08 

2 

-8 

-2.5 

-0.8 

1  .7 

58.7 

60. 1 

7.0 

235.0 

228 . 0 

12.0 

230.6 

218.6 

200 

1  99 

1  1  9 

2 

-8 

-3  .  1 

-1.4 

1  .  7 

58  .  1 

60.7 

8.0 

236.0 

228 . 0 

13.1 

231.6 

218.6 

200 

1  99 

125 

2 

-7 

-3 . 3 

-2.2 

1  .  1 

57 . 3 

61.5 

9.5 

237 . 5 

228 . 0 

14.6 

233  .  1 

2  18.5 

200 

1  99 

143 

1 

-2 

-3.8 

-2.4 

1  .  4 

57.0 

61  .7 

10.0 

238.0 

228 . 0 

15.1 

233.6 

218.5 

20  1 

201 

I  6  I 

1 

3 

-3.9 

-2 . 3 

1  .  6 

57.0 

61.7 

10.0 

238.0 

228 . 0 

15.1 

233.6 

2  18.5 

201 

201 

1  79 

1 

6 

-3  .  5 

-0 . 7 

2.8 

58.6 

60.  1 

7 . 0 

235.0 

228 . 0 

12.2 

230.6 

218.4 

200 

201 

1  97 

1 

5 

-2 . 8 

-0  .  1 

2  :'l 

59.2 

59.5 

5.9 

233.9 

228 . 0 

11.1 

229.5 

2  18.4 

20  1 

201 

2  1  5 

1 

0 

-2  .  5 

C.8 

3  .  3 

60 . 2 

58 . 5 

4.0 

232 . 0 

228.0 

9  .  3 

227.6 

2  18.3 

2  02 

202 

233 

1 

-4 

-2 . 6 

0.0 

2 . 6 

59.4 

59.2 

5 . 4 

2  3  3.4 

228 . 0 

10.7 

2  2  9.  n 

2  18.3 

2  0  2 

20  1 

239 

1 

-4 

-2.6 

-0.3 

2  .  3 

59.0 

59.6 

6.0 

2  34.0 

2  2  8.0 

n  .  3 

229.6 

2  18.3 

2  0  2 

201 

avg  roll  adj=2.2  stdv=0.6  cnt=41 


pij « u  rii ■  u  II  I 
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Ai'PEMDIX  5.  Geographic  Calculations 
" • 1  qeomim 

-'lA.GE;  qeomim  NN  (.AIM  opO  ihH  caA) 
r XAM  geomim  0 1 
IN'FCT  \AIM\f  ilNN.d 
ITFUT  \AIM\geoNN.d 

ipplys  geom  corrn  for  pit  yaw  roll  lat  Ion  alt 
o-nerates  xi  =  k.m  along  orbit, yi  =  km  perp  orbit  and  pathlen 
p '>  1  (debug  test  program  with  interactive  p  y  r  om ) 

:  •  1  H  is  height  of  lorizon  (150  for  ionosph  horizon) 

■".A  A  reduce  scan  angle  by  A  for  high  resolution  reduced  area  (not 
implemented  ) 

generates  a  41  x  41  matrix  of  geographic  km  from 
I'cbi*:  nadir  for  each  image  line  and  column 

1  r m  1  n a  1  output  from  qeomim 

■  F : BIN \GEOMIM . EXE  06 
i  ■  '-ns  - 

i'"  -  '  clipan  =  0  horht  =  0 
■  .  :  :  ri  i  1 1  c  n  s 

..-.a  degrees  from  nadir  to  point 

km,  from  nadir 

k:n  from  first  line  nadir 

:  t  •.■ertical=  l  .0 


•r.  • 

I 

i 

I 


3  0  0  = 


"  3  q  - 


path 


'  ; :  r  inq  = 

'  ■  g  ’  r  1  X  ^  2 

col 

;  r.  q 


pd^  h 


1  "1 1  r  a  r;  g  ■ 
d  t  r  :  X  col 
c  1  q  “  col 
nq 


=  0  i 

mage  lin=  0  p= 

1  y  = 

1 

r  = 

-2  alt=1019 

kmS=  0 

hion= 1 50 

2  .  2 

col=123.6  kmS= 

-12.0 

0 

5  10  15 

20 

25 

30 

35 

40 

1  4 

40  67  93 

1  20 

1  48 

1  76 

205 

233 

-57 

-42  -28  -14 

0 

1  5 

30 

46 

61 

2  005 

-924  -520  -253 

-32 

200 

473 

884 

2005 

-39 

-23  -  1  8  -  1  5 

-12 

-  1  0 

-7 

-4 

3 

2  1 

59  81  94 

99 

96 

83 

61 

21 

20  image  lin= l l 9  p= 

2  y^ 

-8 

r  = 

-1  alt=1027 

kmS=2 1 42 

hion= 1 50 

1  .  5 

col  =  I  22 . 5  kmS  = 

2  116.0 

0 

5  10  15 

20 

25 

30 

35 

40 

1  3 

39  66  93 

120 

147 

175 

203 

231 

-58 

-43  -29  -14 

0 

1  5 

30 

45 

60 

2003 

-917  -512  -245 

-22 

207 

478 

886 

1  994 

2  385 

2239  2185  2149 

2119  2088 

205  1 

1  993 

1  826 

2  1 

59  81  95 

99 

96 

83 

6  1 

2  1 

40  image  lin=238  p= 

1  v  = 

j 

) 

r  = 

-n  al 

1  -  1  P  3  1 

kmS=4284 

h  1  '^n  =  1  5  0 

0  .  3 

col = 1 20 . 6  kmS= 

4263-5 

0 

5  10  15 

20 

25 

3  0 

35 

40 

1  1 

38  65  92 

120 

147 

1  74 

202 

229 

-59 

-44  -29  -14 

0 

1  4 

29 

44 

59 

-202  1 

-916  -505  -233 

-5 

224 

497 

908 

2017 

20 


I."  »  L  (  1.  ■  L'^  L  II  ^ 


i 


“’  ■  w  .■ 


p'j«'  Ml  wm\^  IP 


kmS 

44  1  2 

4331  4301 

4281  4264 

4246 

4225 

4  1  92 

4097 

1  0  0  path 

2  1 

60  82 

95  99 

95 

82 

60 

21 

5.2  so le 1 

USAGE:  solel 

NN 

(AIM  ) 

EXAM:  solel 

01 

INPUTS  \AIR\hdrNN.d  \AIR\geoNN.d 
OUTPUT  \AIR\selNN.d  41x41  sin_sol_el 

Terminal  output 

C:\BOB\BIN\SOLEL.EXE  06 

rolad=  2.20  rot=0 

rho0=  91.5  sg=  0.1  cnt=  30 

lam0=  327.2  sg=  0.6  cnt=  31 

Center  Lat=  69.4  Lon=  326.4 

lalox  lon=145.5  diff=180.9  correcting  by  180 

Day=  47  LMT=  9.86  Solar_noon  lat=  -13.2  lon=  358.5 

def initions 

xi  km  along  orbit  from  nadir  first  line 

yi  km  perp  to  orbit  from  nadir  current  line 

sinse  100  *  sin(solel) 


matrix  samples 
il=  0 


ik 

0 

1 

5 

1  0 

1  5 

20 

25 

30 

35 

39 

40 

sc  ang- 

5  7.3- 

54.4- 

43.0- 

28.6- 

14.3 

0.0 

15.4 

30.8 

46.3 

58.6 

61  .7 

im  CO  1 

1  4 

1  9 

40 

67 

93 

120 

1  48 

1  76 

205 

227 

233 

X 1  /  1  0 

-3 

-3 

-2 

- 1 

-1 

-1 

- 1 

0 

0 

0 

0 

yi/'  1  0 

-200 

-16  1 

-92 

-52 

-25 

-3 

20 

47 

88 

1  59 

200 

lat 

72 

75 

81 

85 

87 

88 

88 

86 

82 

76 

73 

Ion 

24  1 

242 

246 

252 

265 

300 

1  0 

36 

46 

51 

52 

sinse 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

il  =  20 

sc  ang- 

58.1- 

55 . 2- 

43.6- 

29.0- 

14.5 

0 . 0 

15.2 

30.3 

45.5 

57 . 6 

60.7 

im  CO  1 

1  3 

1  8 

39 

66 

93 

120 

147 

175 

203 

225 

231 

X  i  /  1  0 

238 

233 

223 

2  1  8 

2  1  4 

21  1 

208 

205 

1  99 

189 

1  82 

y  i  /  1  0 

-200 

-161 

-91 

-51 

-24 

-2 

20 

47 

88 

1  58 

199 

lat 

61 

63 

67 

68 

69 

69 

70 

70 

69 

67 

66 

Ion 

288 

293 

305 

314 

320 

325 

33  1 

338 

349 

4 

1  2 

sinse 

0 

0 

2 

4 

6 

7 

9 

1  0 

1  3 

1  6 

1  8 

i  1  =  40 

sc  ang- 

59.1- 

56.1- 

44  .  3- 

29.5- 

14.8 

0 . 0 

14.9 

29.8 

44 . 7 

56.6 

59 . 6 

im  col 

1  I 

1  6 

38 

65 

92 

1  20 

147 

1  74 

202 

224 

229 

X  :  '1  'J 

4  4  1 

4  38 

433 

430 

428 

426 

424 

422 

4  1  9 

4  1  3 

409 

"  ^  III 

J  ^ 

-  2  /'  2 

-16  1 

-9  1 

-50 

-2  3 

0 

"> 

40 

90 

161 

20  1 

1  cl  ^ 

46 

47 

49 

49 

50 

5  0 

^11 

(1 

49 

49 

1  n 

3  0  1 

3  06 

3  1  4 

320 

32  3 

i2  7 

3  3  " 

■  cl 

^  3  9 

34  9 

3  54 

c:  1  p;  t-c. 

:.o 

'9  v 

29 

3  2 

34 

35 

37 

3? 

4  1 

45 

47 

S'';'l3r  ele-;  at  image  ctr  =  4.0  sin  =  0.070 


21 


-j»  r^r 


K 

K 

t. 


5  •  \  yj^  yi' 

warpim  NNDp  (AIM  opN  ) 

;eads  selNN.D  and  geoNN.d 
rates  ground  map 

by  applying  geom  corrn  for  pit  yaw  roll  lat  Ion  alt 
NOTE  expects  filtered  image  filNN-iD  in  qui 
opN  n=0  200*sinel  and  200/path  to  qu4 

then  move  q4  to  selNN-i  for  geomlsq 
n=^  I  warps  image  to  q2 
put  rawNND  into  q3 
histogs  and  computes  amplif 
reports  wvln  Rayl  amp  into  \aim\ampNNDp . d 
then  move  q2  to  warNNDp.i  for  geomlsq 
n=2  both  warps 

EXAM  warpim  01  for  solel  image 

EXAM  warpim  012a  opi  for  war012a.i  scene 

Terminal  output 

i' :  \BOB\BIN\WARPIM.EXE  06 

100  *  sin  solar  elev  at  60  line  spacing  0  0  7  26  47 
kmscl  map  si2e 

xysc  scaling  to  240  line  x  240  col  image 

Nctr  km  from  first  line  nadir 

Fctr  km  from  center  line  nadir 

kmscl  5000  xysc  0.048  Nctr  2119  Ectr  -22 


APPENDIX  6.  Geomlsq 


USAGE _ geomlsq  NND  (AIM  hiH  loB  plP  ) 

INPUTS  geoNN.d  in  qu2  selNN.i  in  qu4 
OUTPUT  \AIR\ IsqNND.d 

plP  limit  analysis  to  inv  path  pet  >  P  (def=  25  ==  path  4) 
loB  dont  use  inputs  less  than  B  (def=4) 
hiH  H=l  report  histogrammed  occurrences  of  fit 
EXAM  geomlsq  042 

Terminal  output 

options- 

geomlsq  for  AIR  062  hil 
histoq  of  input 

lut  inverse  path  mn=46  mx=198  (200  =  path  1.0) 
lut  sin(solel)  mn=  0  mx=  93  (200  =  sin ( so  lei )  =  1  ) 

percent  sunlit=63 

reported  as 

p_mn  46  p_mx  198  s_mn  1  s_mx  93  pcsl  63 

pet 

tot  med  99pctl  av  sg  sg/av 

nsIO  his  tot  med  i99  3575  5  196  av  sg  pesg/av  21.1  42.8  202 

slIO  his  tot  med  i99  6160  47  201  av  sg  pesg/av  62.8  42.4  67 

stored  ent  path  sel  for  sampling 
total  getvals  =  4500 

fit  to  log  1 0 ( ent )  =C0  +  Cp  *  log  1 0 ( 1 . /path )  +  Cs  *  log  1 0 ( 200 . *sel  ) 

sdvO  std  dvn  inputs 

sdvf  std  dvn  after  fit 

pevr  perct  variance  remaining 

COi  CO  converted  to  10**C0  and  normalized  to  path=1  sel=0.01 
pcsd  std  dvn  after  fit  converted  from  loglO  to  +  pet  of  value 
pts  points  fitting  within  +  limit 


first  fit 
fc  si  3.02 
screened  fit 


firs+'.  fit 
Is  s 1  2.97 

screened  fit 
1  si  2.77 


sdvO 

:p=  -1 

sdvf 
Cs  =  0  . 

pevr 

5 

COI 

Cp 

Cs 

pcsd 

pts 

0.17 

0.15 

73 

7.39 

-1.00 

0.50 

4  1 

1  675 

0.23 

t 

0 .09 

1  4 

6.55 

-1.00 

0  .  50 

22 

2449 

0.17 

0.13 

53 

14.37 

1 

o 

00 

0.27 

34 

1  675 

0  .22 

0 . 09 

1  5 

P  .  38 

-0.84 

0.41 

22 

2466 

h,  I  ‘ ;  t  ■ '  cj  r  a  m  o  f  s  a  [ii  p  1  e  s 
^;unllt:  samples -=822 

\'-:n  br  below  =  log  (  br  )  -  0  .  5*  log  (  se  1  ) 
no_s-'^'  samples  =  822 

h.'s^-oq  cell  above  9  (A=10-19  B=20-29 


log  p i n  V 
path 

1  g  b£  ^ 
.to  4 


AE. 

.  3  2ABCD. 
n28BACBBE. 
1 57BBBAAAAAA. 
.  2  1 1 2AABB93765757A. 

I  112  1BAAA64A56234167. 
11  12137553AA7964732714385. 

II  1 5585547A79651 1 223346A. 


’  .  O'  0 

1  G  . 

.123  1 

3365334856325 

.22  3  1212 

1  1  . 

1  221 53545463 

3341  41 

.2  211  2 

1  3  . 

2263651747.22 

141  1 

.12  11 

'  .  ■  5 

'  4  . 

2343581 1312 

2 

113111  1  1 

■ 

1  6  . 

445533 

1 

1 

1  1  1 

.  _  4 

1  P  . 

52121 

.  1 

1  1 

.11  111 

'  .  3  0 

2  0  . 

2  1 

1  . 

1 

1  111 

'  .  35 

22  . 

1 

.  1 

1  1 

'  .  4  ^ 

25  . 

• 

1  1 

1 

1.45 

28  . 

.  1 

12  1 

.  1 

• 

r  ■■  c  3  'a  n 

fit 

( Ig  ICO 

sdvO 

sdvf 

pcvr 

C01 

Cp 

Cs 

first 

fit 

Is  ns 

1.01 

0.16 

0.16 

1  00 

6.08 

-0.10 

0.00 

screened  fit 

Is  ns 

1.03 

0.13 

0.13 

99 

5.80 

-0.12 

0.00 

Cs  pcsd 


h i stog rammed  samples 
no__sun  samples=l89 
rio__s=l  samples=l89 

hipt.og  cell  above  9  (A=10-19  B=20-29  ..) 


*0'g  pinv 
path 
i  r  b  r 

' .  6  0'  4 


163895966324334A. 
.  212545553216424687. 

1 .  1234  343432343A. 

23  113131  4251312377. 

.  22  332  2  11  31 2427 . 

2  21 1 32  2.  1  13155. 

2131321.111 324  . 


iag  1  ng 


with  aur  thr=47  slsgrat=143  dksgrat=171 


■  ij-  ^  tograp:  after  reimaging 

bnr  his  tot  med  i99  19931  42  143  av  sg  pcsg/av  5  1.8 
ns>br  his  tot  med  i99  12529  4  31  av  sg  pcsg/av  8.5 
o’br  his  tot  med  i99  4396  117  224  av  sg  pcsg/av  124.1 


28.2  54 

8.5  99 


his  tot  med  i99  4396  117  224  av  sg  pcsg/av  124.1  53.7  43 
his  tot  med  i99  19929  101  199  av  sg  pcsg/av  107.3  24.2  22 
his  tot  med  i99  12531  99  239  av  sg  pcsg/av  123.0  82.1  66 


24 


opening  \AIR\ lsq062 .d  for  output 


AIR 

Site 

Day  GMTC 

LMT 

Linp 

Late 

LonC 

Dir 

062 

SON 

047  12. 

1  0 

9.86 

229 

69.4 

326.4 

S 

Rad  j 

Wvl 

Rayl  mxse 

Ampk 

Araw 

Af  il 

pcsl 

pcau 

2 . 2 

1356 

30 

48 

1  .  000 

7.82 

46.7 

63 

1  2 

source 

pts 

avg 

sig 

med 

PC_ 

sig/avg 

s  Ibr 

1  993  1 

5  1.8 

28.2 

42 

54 

nsbr 

12529 

8 . 5 

8 . 5 

4 

1  00 

aubr 

4396 

124.1 

53.7 

1  1  7 

43 

sift 

1  9929 

107.3 

24 . 2 

101 

23 

nsf  t 

12531 

123.0 

82.  1 

99 

67 

APF'FNDIX  7.  Multiband  Process 


ir. u  1  bnd  help 


USAGE  mulbnd  stS  blB  whW  opX  klK  khK  IIL,  IhL  bgB 
scales  inputs  qul2(3)  to  output  screen  4  coded-color 
cpX  use  in  2-color  inputs 

X=l  shows  color  without  brightness  scaling  X=2  same  with  black  for 
zero  vals 

reports  histograms  from  klK  to  khK  and  IIL  to  IhL 

stS  S=0  no  stretch  (each  assumed  pre-scaled  0-255)  S  >  0  stretch  S  pctl  to 

254*3/100 

S  <  0  histog  equ  with  stretch  cmpc**abs(S) 
bgB  'bug'  special  use  in  2-color  inputs 
B  =  2  reports  input  hists  output  stretch  luts  and  quits 
B  =  1  generates  stretched  quad  1  and  2  onto  quad  3  and  4  count  only 
rrom  blB  (def=0)  to  whW  (def=255) 

EXAM  mulbnd  qu1204  (default  2-color  1  and  2  to  output  4)  EXAM  mulbnd  qu2143 
(3-color  2  1  4  to  output  3) 

Terminal  output 

'c \BOB\BIN\MULBND.EXE 
options- 

:’’>ulbnd  fi=  1  2  0  fo  =  4  op  =  0  st  =  0  bg  =  0 
histogram  report 

band=i  samples  0=10362  1=15436  255=160  histog  other=28800 
band=2  samples  0=10362  1=16273  255=0  histog  other=28800 
histog  of  bandl  vs  band2  values  mx=122  excluding  O’s 

band2  0  15  30  45  60  75  90105120135150165180195210 

h)  n  'i '' 


0 

99 

0 

2 

99 

55 

37 

38 

28 

19 

21 

25 

1  7 

7 

3 

0 

1  5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

45 

99 

0 

0 

1  8 

4 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

60 

99 

0 

4 

66 

40 

4 

1 

0 

3 

6 

0 

1 

0 

0 

0 

7  5 

99 

0 

0 

53 

32 

1  2 

0 

1 

0 

1 

2 

0 

0 

0 

0 

9  0 

75 

0 

0 

1  9 

1  8 

25 

0 

4 

1 

3 

4 

1 

1 

0 

0 

;  0  5 

56 

0 

1 

22 

1  5 

1  5 

1  2 

6 

0 

1 

1 

0 

0 

0 

0 

1  20 

55 

0 

3 

1  1 

1  7 

1  9 

6 

6 

4 

9 

2 

4 

0 

0 

0 

1  35 

6  3 

0 

0 

1  5 

1  4 

1  5 

7 

7 

5 

6 

1  0 

1 

0 

0 

0 

1  5  0 

47 

0 

0 

1  2 

1  0 

4 

4 

1  1 

6 

6 

16 

3 

0 

0 

0 

'  ^■>  5 

1  7 

0 

0 

9 

1  6 

22 

2 

9 

1  4 

1  7 

9 

9 

1 

0 

0 

1  80 

7 

0 

0 

9 

1  3 

4 

7 

16 

22 

1  1 

1  7 

1  0 

5 

0 

0 

195 

0 

0 

0 

1  3 

4 

3 

6 

9 

9 

1  3 

3  1 

26 

6 

0 

0 

/  ’ 

n 

0 

0 

1 

4 

6 

2 

1  9 

33 

34 

94 

99 

50 

26 

0 

Mc'te  to  equalize  signal  response: 


histogram  equalize 


band  i 


cmpc 

0 

1 

2 

1  0 

50 

80 

90 

95-100 

1  va  1 

0 

64 . 9 

65 . 2 

69.0 

1  35 

222 

248 

255 

1  ut 

0 

3  .  5 

4 . 7 

2  1.1 

127 

204 

229 

255 

band  2 


-‘•'A',*.’ '.'.•-Vv’ -v  AIiMiAlAj.  •> 


cmpc 

0 

1  2 

1  0 

50 

80 

90 

95 

98 

99 

99.5 

99.8 

1  00 

i  val 

0 

51.0  53.1 

59.0 

1  22 

182 

1  97 

206 

2  1  9 

222 

224 

234 

237 

lut 

0 

1.0  4.4 

23.1 

1  27 

204 

228 

242 

250 

252 

254 

255 

255 

mulbnd  of  histoqram-eaualized  bands 


Note  improved  distribution  of  data 


band=l  samples  0=25854  1=0  255=160  histog  other=28800 
band=2  samples  0=26643  1=9  255=23  histog  other=28800 
histog  of  bandl  vs  band2  counts  mx=65  excluding  O's 

band2=0  15  30  45  60  75  90105120135150165180195210 


band  l 


0 

99 

86 

69 

50 

59 

38 

84 

61 

1  5 

99 

6 

2  1 

1  6 

6 

0 

0 

1 

3  0 

99 

7 

4 

7 

9 

3 

3 

0 

45 

99 

1  3 

1  5 

1  9 

1  2 

3 

0 

0 

60 

99 

1  3 

6 

1  5 

1  8 

9 

1 

1 

75 

99 

2  1 

35 

1  9 

1  3 

1  2 

3 

1  2 

90 

99 

1  5 

27 

1  6 

2  1 

43 

1  5 

1 

105 

99 

7 

1 

1  6 

1  8 

24 

4  1 

7 

1  20 

99 

4 

9 

1  5 

26 

29 

26 

7 

1  35 

96 

4 

1  8 

1  0 

1  0 

4 

1  2 

27 

1  50 

58 

6 

1  0 

4 

32 

2  1 

24 

1  6 

165 

1  9 

6 

1  2 

24 

7 

9 

7 

36 

’80 

3 

1  2 

6 

1 

0 

6 

6 

33 

195 

1  3 

0 

0 

0 

6 

0 

6 

6 

2  1  0 

0 

0 

0 

0 

0 

0 

0 

0 

1  0 
0  0 
0  4 

1  0 
2  3 

1  3 

7  7 


39  23  33 
1  0  0 


0  0 
4  0  1 

0  7  0 

4  1  0 

0  4  10 

4  3  1 

10  13  6 

10  16  19 
38  12  10 
35  1529 
35  29  32 
26  46  44 
23  7  35 


4  1  0 
3  1 


16  24  21 
0  0  0 
0  0  0 
0  0  0 
4  3  0 

0  0  0 

3  0  3 

4  3  0 

7  1  0 

4  3  0 

26  6  6 
16  12  13 
43  46  16 
66  39  99 
52  50  63 


A ;  ;  t-lN'D I X  8.  Batch  Processes  MULP  and  MULD 
■  of  MULP.BAT 

'  ' " ' ' ji i ^  1  2 

■  :"i;lp  prepares  2-band  colored  images  from  raw\air  images 
;  ■’  "■  ild  demonstrates  using  images  prep  by  mulp  for  fast  demo 

processes  via  unpak  filtraw  warpim  geomlsq  mulbnd 
rerun  can  bypass  parts  with  'goto  aa  etc' 
r-^m  final  images  puts  landmass  squishes  to  square  and  titles 
rem  CR  if  raw%l%2  raw%1%3  OK  else  CNTRL  C 

bs^ch  command  purpose 


It  not  exist  \air \raw% 1 %2 . i  goto  erri 
if  not  exist  \air \raw% 1 %3 . i  goto  erri 

modlut  cw 

rem  geomim  solel  warpim  opO  once  only 
rem  serves  all  bands 


test  if  present 

fl 

set  display  to  4-wedge  color 


raw  fit 

%1 

prepare  fil%1.d  orbit  and  hdr%1.d 

geomim 

%  1 

prepare  40x40  matrix  km  geoNN.d 

^  ■  -'lei  %  1 

prepare  40x40  matrix  solel  selNN.d 

w  1  p  i  m. 

%  1 

prepare  image  lut  of  path  and  sel 

moveim 

45 

\air \sel% 1 . i 

store 

rem  now 

proc  %  1 %2 

n-'  im 

52 

\air\raw% l %2 . i 

xfer  raw%1%2  image  to  q2 

.!'.pak  2 

1  % 

1  %2 

unpack  packed  values 

'  .'e im 

1  5 

\air \unp% 1 %2  .  i 

store 

i:  12  - 

_  M 

1 "  filtraw  n7 

filter  raw  image  (det  1  wid  =  7  ) 

if  not 

%2 

==  "1"  filtraw 

n3 

filter  raw  image  (det  2,3  wid=3) 

r'O'.'e  im 

45 

\air\fil%i%2.i 

store 

IT!  e  i  m 

5  3 

\air\unp% l %2 . i 

move  unp  to  q3 

mo  ■/e  im 

5  1 

\air\fil%l%2.i 

move  fil  to  ql 

warpim 

%  1  %  2  op  1 

warp  to  q3  and  compute  amplif 

mo ve  im 

25 

\air \war% l %2  .  i 

store 

r  ;  e  1  m 

52 

\dir\war%i%2.i 

move  warped  to  q2 

•  1  ni 

54 

\air \sel%  l  .  i 

move  sel  lut  image  to  q4 

%1 

%2 

Isq  fit  to  path  and  sel 

pic'/eirn 

1  5 

\air \agl% 1 %2  .  i 

store  airglow 

rno  im 

3  5 

\ a ir \ aur  %  1 %2  .  i 

store  aurora 

--'m  now 

proc  %  1  %  3 

rr^  ■ ; '/  e  i  rn 

52 

\ a ir \ raw% i % 3  .  i 

repeat  for  NND  D=%3 

■o.pak  2  1  % 

1 

mo  *.'e  :  ni 

1  5 

\air\unp%l%3.i 

_  »» 

1”  filtraw 

n7 

\  f  n  o  i 

b  <L. 

==  "1"  filtraw 

n  3 

i  m 

4  ^ 

'air\fil%l%3.i 

T  "'  rT’ 

53 

\ a  1 r \unp% 1  %  3 . i 

'  pi 

1  1 

'  a  i  r  \  f  i  1  %  1  o  3  .  i 

w  j  r  ['■  i  m 

i  I  1 

■  opl 

P'1 

2  5 

'  .  a  i  r  \  w  a  r  %  1  ?■  3  .  i 

r  '/fM  m 

5  2 

\ a i r \ war^ I  %  3  .  ] 

P  '  1  rn 

54 

\ajr\sel'^l  .i 

J-  '  ’n !  eq 

1 

%3 

28 


•»TT  J1.J1 


'  V  -  V  ■■  V^THTT"^  'V  VI  V- 


moveim  15  \air \agl% 1 %3 . i 
moveim  35  \air\aur%l%3.i 
moveim  53  \air\aur%1%2.i 
moveim  54  \air\aur%l%3.i 
h 1 seq  31  hi  b2 
hiseq  42  hi  b2 
mu Ibnd 

if  exist  \air\lm%i.l  goto 
landmass  q3  a%l 
moveim  35  \air\lm%l.i 
:  Imx 

moveim  53  \air\lm%l.i 

overlim  31 

cverlim  32 

o e  r  1  i  m  3  4 

moveim  1 3 

rot im  31  0  1.2  1 

moveim  23 

r ot im  32  0  1.2  1 

moveim  43 

rotim  34  0  1.2  1 

texthdr  % l %2  1 

texthdr  %1%2  2 

texthdr  %I%2  4 

texthdr  %1%2  1  1 

texthdr  % l % 3  2  3 

text hdr  %  1 %2  4  2 

texthdr  %  1  % 3  4  4 


prepare  for  multi-bnd  aurora 


hiseq  aurora 
prepare  coded  image 

prepare  landmass  overlay 
store 

recover 

overlay  landmass  on  aur  qd2 
"  qd3 

"  qd4 

prepare  'squished' 


texthdr  %1%2  2 

texthdr  %I%2  4 

texthdr  %1%2  1  1 

texthdr  % l % 3  2  3 

text  hdr  %  1 %2  4  2 

texthdr  %  1  % 3  4  4 

moveim  15  \air \dau% 1 %2 . i 

moveim  25  \air \dau% 1 %3 . i 

moveim  45  \air \dau% 1 %2%3 . i 

rem  now  display  warped  airglow+aurora 

moveim  53  \air \war% 1 %2 . i 

moveim  54  \air \war% 1 %3 . i 

hiseq  3 l  hi  b2 

hiseq  42  hi  b2 

mu Ibnd 

if  exist  \air\lm%i.i  goto  Imy 
landmass  q3  a%l 
moveim  35  \air\lm%1.i 
;  Imy 

moveim  53  \air\lm%l.i 
overlim  31 
overlim  32 
o'/erlim  3  4 
exthdr  %  1  %2  1 
♦■‘=-x''hdr  ^i"2  2 
t^x^  hdr  3  1  i  2  4 
‘  c.x''  hdr  3  1  i  2  1  1 

X  h  ':i  r  ^  i  ;•(  2 
'■-•x'’hdr  '  'i>2  4  2 

’  X  ^  h  d  r  %  1  %  3  4  4 
V'  '/eLm  15  a  1  r  '  dwr  %  1  % 2  .  1 
m'  1  m  2  5  ;  a  i  r  \  dw  r  %  1  ^  3  .  i 


headers  on  quads 


\air\lm%i.i  goto  Imy 
q 3  a%  1 
\air\lm%1  .i 

\air\lm%l  .i 


store 


recover  warped  filtered  images 

hiseq 

mu Ibnd 

prepare  landmass  overlay 
store 

get  landmass 
overlay 


store 


ei[!i  45  \air\dwr%l%2%3.i 
!  his__eq  and  mulbnd  scene  NN  D  D 
trom  raw  and  filtered  images 
5  3  \  a i r  \  raw%  1  % 2  .  i 
3  ’  hi 
2  s  2 


'4 


‘.‘o  1  rr'. 

■  .7  4 

55 

• 

.■  \ 

7;'.:bnd 

-  -x^  hdr 

*  y  *-  ^7  r- 

‘  x"  hdr 

*  >;  *■  ■:  m 

H 

■  t 

m  ■ :  e  i  in 

•V 

1 .  • 

rem.  now 

m  '  e  i  m 

iii.seq  4 

.H' 

■  ;  rn  i  4 

i'  criq  q'  hiseq  q2  ampl 
"x*:  im  q'j3  1134  fml  orig  raw%l%2 
c<tim  qu4  1134  fml  hiseq 

im  qu2  1134  fml  stretched 
c<-hdr  ^i32  3 
■x^hdr  3  1 

t  c  r  e 

25  \  a  ir \drs%  1  %2  .  i 
. ■.■eim  45  \air\drh%l%2.i 
;''.eim  35  \air\drl%l%2.i 


get  rawNND  D= 1 

histog  equalize 

stretch  2-998  permil  0-255 


overlay  texts 


store 


repeat  for  rawNND  D=%3 


mulbnd 


1  %2  4 


store 

repeat  for  filtered  images 


h::-''=‘q  4  2  hi 
m  ■  j  1  b  n  d 

♦  ^-xT  irn  qu4  1134  fml  fil  hiseq  2-band 
:n-.  /eim  45  \air\df2%1%2%3.i 

qo"o  end 

:  >*  r  1 

•  raw'll  %2  or  %1%3  not  found 
;  -c. i 
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MULD.bat 

rem  USAGE  muld  NN  D  D 

rem  mulp  prepares  2-band  colored  images  from  raw  \air  images 

rem  muld  demos  using  images  prep  by  mulbp  for  fast  demo 

rem  processes  via  unpak  filtraw  warpim  geomlsq  mulbnd 

rem  rerun  can  bypass  parts  with  'goto  aa  etc' 

rem  final  images  puts  landmass  squishes  to  square  and  titles 

rem  CR  if  raw%l%2  raw%1%3  OK  else  CNTRL  C 

rem  pausex 

20omim 

modlut  cw 

rnoveim  5  1  \ air  \  r aw%  1  %2  .  i 

h 1 seq  12  s  1  0 

texfim  qul  sz2  fml  raw%l  det%2 
textim  qul  1134  fml  orig 
textim  qu2  1134  fml  stretched 
rnoveim  53  \air\fil%l%2-i 
textim  qu3  1134  fml  filtered 
rnoveim  54  \air \war% 1 %2 . i 
textim  qu4  1134  fml  geom  warped 
pausex 

rnoveim  51  \air \agl% 1 %2 . i 
textim  qul  S22  fml  raw%l  det%2 
textim  qul  1134  fml  Isq  airglow 
rnoveim  52  \air \aur% 1 %2 . i 
textim  qu2  1134  fml  Isq  aurora 
pausex 

rnoveim  51  \air \raw% 1 %3 . i 
h i seq  12  s 1 0 

textim  qul  sz2  fml  raw%1  det%3 
textim  qul  1134  fml  orig 
textim  qu2  1134  fml  stretched 
rnoveim  53  \air\fil%1%3.i 
textim  qu3  1134  fml  filtered 
rnoveim  54  \air \war% 1 %3 . i 
textim  qu4  1134  fml  geom  warped 
pausex 
:qq 

rnoveim  51  \air \agl% 1 %3 . i 
textim  qul  sz2  fml  raw%1  det%3 
textim  qul  1134  fml  Isq  airglow 
rnoveim  52  \ a ir \ aur %  1  % 3 . i 
rextim  qu2  1134  fml  Isq  aurora 
pausex 

rnoveim  5'  \ a i r \dwr %  1  % 2 . i 

‘v'Xtim  qul  1134  fml  aur  +  airg  det  %2 

nv'  1 IV  5  r,  \  a  1  r  \  d  w  r  1  %  3  •  i 

f im  qui:  1134  fml  aur  +  airg  '■l‘-''t  "-3 

I  •v’l  rv^-w  display  hiseq  aurora 

-V  ■;  i  II  5  j  \  a  1  r  d  a  u  %  1  %  2  .  1 

qul  1134  fml  Isq  aurora  det  %2 
m  ► '  1  m  5  4  \  a  i  r  \  d  a  u  %  1  %  3  •  i 


textim  qu4  1134  fml  Isq  aurora  det  %3 

pausex 

nv'd  l\.it  c2 

n’o-.’eim  51  \  a i r  \dwr  %  1  %2  % 3  .  i 

taxtim  qul  sz2  fml  2- band  geom 
moveim  52  \air \dau% 1 %2%3 . i 
r.extim  qu2  sz2  fml  2-band  aurora 
r  a  u  s  o  x 
z  com  im. 
modlut  cw 
:  read  from  disk 
r'oveim  51  \ air \dr  1  %  1  % 2  .  i 

mcveim  52  \air \drh% 1 %2  .  i 
moveim  53  \air \drs% 1 %2  .  i 
moveim  54  \air \ f i 1% 1 %2 . i 
textim  qu4  1134  fml  filtered 
pausex 

rem  black-white 
modlut  cO 
pausex 
modlut  c2 

moveim  51  \ a ir \dr 2  %  1 %2 %  3  .  i 

moveim  52  \air\df2%l%2%3.i 

move_m  53  \air \dwr% 1 %2%3  .  i 

textim  qu3  1134  fml  geom  aurora+airglow 

moveim  54  \air \dau%  1  %2%3  .  i 

textim  qu4  1134  fml  geom  aurora  only 


APPENDIX  9.  Photographs 


A  series  of  color  photographs  have  been  taken  from  the 
monitor  to  document  processing  steps.  The  color  images  cannot  be 
reproduced  by  an  office  photocopier  but  4  sets  have  been  prepared 
for  4  copies  of  this  report. 
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0-63  ->  red  64-127  ->  yellow  128-191  ->  green  192-255  ->  blue 
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